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ABSTRACT 

Context. Nearly 100 new strong lens candidates have been discovered in the COSMOS field. Among these, 20 lens candidates with 
0.34< ziens <1.13, feature multiple images of background sources. 

Aims. Using the multi- wavelength coverage of the field and its spectroscopic follow-up, we characterize the evolution with redshift 
of the environment and of the dark-matter (DM) fraction of the lens galaxies. 

Methods. We present spectroscopic and new photometric redshifts of the strong lens candidates. The lens environment is characterized 
in the following way: we account for the projected 10 closest galaxies around each lens and for galaxies with a projected distance less 
than 1 Mpc at the lens galaxy redshift. In both cases, we perform similar measurements on a control sample of "twin" non-lens early 
type galaxies (ETGs). In addition, we identify group members and field galaxies in the X-ray and optical catalogs of galaxy groups 
and clusters. From those catalogs, we measure the external shear contribution of the groups/clusters surrounding the lens galaxies. 
The systems are then modeled using a SIE for the lens galaxies plus the external shear due to the groups/clusters. 
Results. We observe that the average stellar mass of lens galaxies increases with redshift. In addition, we measure that the environment 
of lens galaxies is compatible with that of the twins over the whole redshift range tested here. During the lens modeling, we notice 
that, when let free, the external shear points in a direction which is the mean direction of the external shear due to groups/clusters 
and of the closest galaxy to the lens. We also notice that the DM fraction of the lens galaxies measured within the Einstein radius 
significantly decreases as the redshift increases. 

Conclusions. Given these, we conclude that, while the environment of lens galaxies is compatible with that of non-lens ETGS over a 
wide range of redshifts, their mass properties evolves significantly with redshift: it is still not clear whether this advocates in favor of 
a stronger lensing bias toward massive objects at high redshift or is simply representative of the high proportion of massive and high 
stellar density galaxies at high redshift. 

Key words. Cosmology: Gravitational lensing 



1. Introduction 

In the field of gravitational lensing, studies of strong galaxy- 
galaxy lenses have recently become possible on statistical 
grounds. Indeed, within just a few years, searches in the SDSS 
(SLACS: Bolton et al. 2006, 2008, Allam et al. 2007), in 
COSMOS (Faure et al. 2008, hereafter Paper I; Jackson 2008) 
and in the CFHT-LS surveys (SL2S: Cabanac et al. 2007, 



* Based on observations made with ESO telescopes at Paranal obser- 
vatory under program ID 077.A-0473(A) and under large program ID 
175.A-083911. 



Limousin et al. 2009b) have delivered more than two hundred 
new strong galaxy-galaxy lenses, spanning a wide range in red- 
shift and image angular separation. The reasons for this interest 
are several. First, strong lens galaxies provide measurements of 
the total mass distribution on galaxy scales, bringing additional 
information about the processes of galaxy formation and evolu- 
tion (e.g. Ofek et al. 2003, Chae et al. 2006, Koopmans et al. 
2006, Tortora et al. 2010). Second, from a statistical point of 
view, strong lensing occurrences trace the abundance and con- 
centration of matter in the universe, providing another test for 
cosmological models (e.g. Keeton 2001, Bartelmann et al. 2003). 



1 



C. Faure et al.: Environment and mass properties of COSMOS lens galaxies. 



Studies of the physical and environmental properties of low 
redshift strong lens galaxies (z<0.3) have shown that they are 
bona-fide massive ETGs (Treu et al. 2006, 2009 -T09 hereafter-, 
Gavazzi et al. 2007, Grillo et al. 2009). The strong lens sample 
in the COSMOS field (Paper I), spans a higher redshift range 
(0.33<z<1.13) and previous work on this sample have shown 
that the projected distribution of lens galaxies in COSMOS is 
comparable to that of ETGs, whether in rich environments such 
as cosmic filaments, or in the field (Faure et al. 2009, hereafter 
Paper II). Interestingly, this result means that multiple-mass- 
sheets do not contribute significantly to make more efficient 
lenses, contrary to results found in ray tracing through numer- 
ical simulations (Wambsganss et al. 2005, Hilbert et al. 2007, 
2008). In Paper II we also demonstrate that the presence of large 
scale structures (LSS) in the lens galaxy environment increases 
the angular separation of the lensed images of a source, as no- 
ticed earlier by Oguri et al. (2005). 

In addition, analyses of the mass density profiles of strong 
lens galaxies from hydrodynamical N-body simulations (Dobke 
et al. 2007, Limousin et al. 2009a) and lens samples (Limousin 
et al. 2007, Auger 2008, Natarajan et al. 2009) show that the 
mass distribution in galaxies depends on their environment. For 
example, compared to galaxies centrally placed in their group 
or cluster, some skewing of the total mass density slope appears 
in galaxies located at the edges of a group or cluster. This is 
consistent with tidal stripping of their dark matter halo (T09). 

Therefore, it is important to explore the local environment 
of the COSMOS strong lenses rather than LSS as in Paper II. 
We also want to learn about the lens galaxy mass properties, 
in particular to derive information relative to their DM fraction 
(fz)M, defined as foAf=l-M*/(M*+M£)Af). This is possible as the 
COSMOS strong lenses pertain to one of the best studied region 
of the sky, for which deep multi-wavelength observations have 
been collected (Koekemoer et al. 2007, Taniguchi et al. 2007, 
Capak et al. 2007a, Scoville et al. 2007, Lilly et al. 2007, Ilbert 
et al. 2009). Thus the properties of the COSMOS galaxy popula- 
tion have been studied extensively (e.g. McCracken et al. 2007, 
Capak et al. 2007b, Scarlata et al. 2007a, 2007b, Ilbert et al. 
2010). 

The paper is organized as follows. In § [2l we provide im- 
proved redshifts and stellar masses of the COSMOS lens galax- 
ies. In § [3} we examine and discuss the lens environment. In §|U 
we present results from our strong lens mass models for the sub- 
sample of triple and quadruple image systems. Discussion and 
conclusions appear in § [5j Throughout this paper, we assume a 
WMAP5 ACDM cosmology with D m = 0.258, £2 A = 0.742, 
H = 72 h 12 km s" 1 Mpc" 1 . 



has discovered two additional convincing strong lens candi- 
dates producing multiple images: J095930. 93+023427. 7 and 
J100140. 12+020040.9. The lens galaxies both appear to be 
ETGs. 

Among the 21 multiple image lenses now available, one 
system (COSMOS 5921+0638) - being the only confirmed 
lensed quasar of the sample - has been studied separately 
by Anguita et al. (2009, hereafter Paper III). Another system, 
COSMOS 5737+3424, is a galaxy cluster that is lensing a set of 
background galaxies. In such a case, the lens covers a different 
mass scale than in the other 20 systems, and the lens potential 
is more complex than a single lens galaxy: this target has been 
dropped from the present analysis. 

2.1. Spectroscopic followup 
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Fig. 1. VLT/FORS1 spectra of three strong lens candidates. The 
spectral resolution is 5.52 A per pixel, and has been smoothed 
by a 3 pixels box to improve the quality of the display. Flux is in 
arbitrary units. 



2. The COSMOS sample of strong lenses: new 
redshifts and stellar masses 

By visual inspection of stamp images of 10"x 10" around -9500 
early photometric type galaxies (with redshifts : 0.2 < z < 1.0, 
and absolute magnitude: M v < -20) in the COSMOS field, we 
have discovered 60 strong lens candidates (Paper I). In addi- 
tion, in the same study, 7 lens candidates were found serendip- 
itously across the field. Among the whole sample, 19 systems 
display long curved arcs or multiple images with similar colors 
(based on SUBARU images, by Taniguchi et al. 2007) around 
a bright lens galaxy, and with an image arrangement around the 
lens galaxy which is consistent with the lensing hypothesis, as 
probed by lens modeling. 

Independently, Jackson (2008) inspected the complete set 
of galaxies in the COSMOS field, without discrimination. He 
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- z=0.976±0.001 
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Fig. 2. Keck/DEIMOS spectrum of COSMOS 0018+3845. The 
spectrum is smoothed by 2. 3 A and binned down to 7. 5 A per 
pixel, taking the wavelength dependent noise into account, to 
improve the apparent signal-to-noise. Flux is in arbitrary units. 
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On March 3 rd and April 25-26 th 2006, we have successfully 
obtained spectra for five of the strong lenses in the COSMOS 
field (PI: Faure) using the FORS1 instrument (ESO/VLT) in 
multi-object spectroscopic mode with 19 slits. For each tar- 
get, the central slit was located on the lens candidate with an 
orientation intersecting both the lens galaxy and the brightest 
lensed image. The other 18 slits in the ~7'x7' field around 
the lens, were preferentially positioned on galaxies with colors 
similar to that of the lens galaxy (likely at the same redshift). 
Leaving aside COSMOS 5921+0638 (discussed in Paper III) 
and the lens cluster COSMOS 5737+3424, we display the spec- 
tra of the other three lens galaxies in Fig. [T] Regarding the 
source redshifts, from this dataset we could extract only one 
source redshift, in COSMOS 5921+0638 (Paper III). Positions 
and redshifts of secondary targets around COSMOS 0012+2015, 
COSMOS 0013+2249 and COSMOS 0049+5128 are given in 
Tabled 

In addition, several COSMOS lens candidates have been ob- 
served with the Keck telescopes. Using the Deep Extragalactic 
Imaging Multi-Object Spectrograph (DEIMOS) a spectrum of 
COSMOS 0018+3845 has been obtained on February 11, 2010 
(Fig. 0). The data were collected in 7x1 800s exposures un- 
der photometric conditions with 0.5-0.7" seeing. We used a 1" 
slit together with the 830 line/mm grating tilted to 7860A and 
the OG550 blocking filter. The resulting spectral resolution is 
< 3. 3 A, depending on seeing and object morphology. The ob- 
jects were dithered along the slit by +3" between exposures to 
improve background subtraction, as described in Capak et al. 
(2009). The redshift of the lens has been identified using the 
Call H & K absorption lines and G-band. The source redshift is 
estimated based on the Lyman break which is somewhere in the 
Subaru IA624-band (A effective =6226 A, FWHM=299A), putting 
constraints on the source redshift: 3.9<z s <4.1. From there, we 
find that for z,=3.96+0.02, SilV and CII would align with spec- 
tral absorption features of the lens. This estimation of the source 
redshift is used throughout the paper. 

Moreover, with the Low Resolution Imaging Spectrometer 
(LRIS), Lagattuta et al. (2010) measured the redshifts for seven 
lens galaxies and two sources. 

Finally, from the z-COSMOS catalog (Lilly et al. 2009) we 
have retrieved the redshifts of seven strong lens candidates from 
our list (based on the G-band and Mgl absorption lines). Two of 
them, re-observed with Keck/LRIS, benefit from a second red- 
shift determination from a higher signal-to-noise spectrum. Both 
redshifts are in agreement. 

A summary of the spectroscopic observations available for 
this sample is given in Table [T] while the lenses' redshifts are 
provided in Tabled 

2.2. Improved photometric redshifts 

Recently a new set of photometric redshifts for galaxies in the 
COSMOS field has been derived by Ilbert et al. (2009), using 
the code Le Phare (S. Arnouts & O. Ilberfl) and the photomet- 
ric dataset by Capak et al. (2007a and 2009). Compared to the 
previous set of COSMOS photometric redshifts by Mobasher et 
al. (2007), which relied on 8 broad bands, the main improve- 
ment is the use 30 bands: the 8 broad bands plus 12 intermediate 
bands, 2 narrow bands, 7 bands in the near-infrared and 1 band 
in the ultra-violet. In Fig. [3] we have displayed a comparison 
between the new photometric redshifts and spectroscopic red- 
shifts for the 17 lens galaxies in our sample for which we have 

1 Available at www.oamp.fr/people / arnouts /LE_PHARE.html 



Table 1. Summary of the VLT (FORS1 and VIMOS) and Keck 
(DEIMOS and LRIS) observations: target name, instrument and 
exposure time (in ks). 



Name 


Instr. 


E.T. 






ks 


COSMOS 0012+2015 


FORS1 


1.8 


COSMOS 0013+2249 


FORS1 


1.2 


COSMOS 0049+5128 


FORS1 


1.2 


COSMOS 5921+0638 


FORS1 


1.8 


COSMOS 0050+4901 


VIMOS 


3.6 


COSMOS 0056+1226 


VIMOS 


3.6 


COSMOS 5947+4752 


VIMOS 


3.6 


J095930.93+023427.7 


VIMOS 


3.6 


COSMOS 5857+5949 


VIMOS 


3.6 


COSMOS 0124+5121 


VIMOS 


3.6 


COSMOS 0227+0451 


VIMOS 


3.6 


COSMOS 0018+3845 


DEIMOS 


12.6 


COSMOS 0038+4123 


LRIS 


3.6 


COSMOS 0050+4901 


LRIS 


5.4 


COSMOS 0056+1226 


LRIS 


5.4 


COSMOS 0211 + 1139 


LRIS 


5.4 


COSMOS 0216+2955 


LRIS 


3.6 


COSMOS 0254+1430 


LRIS 


1.2 


J100140. 12+020040.9 


LRIS 


5.4 



Notes. The dispersion is 5.52 A per pixel for FORS 1/1501 and 2.50 A 
per pixel for VIMOS/MR/OS-red. The VIMOS redshifts were obtained 
from the z-COSMOS follow-up of the field (Lilly et al. 2007). The 
Keck/DEIMOS observations were taken with the 830 line/mm grating, 
and using the OG550 blocking filter, leading to a dispersion of 0.47 A 
per pixel. The Keck/LRIS targets were simultaneously observed with 
both blue grating (300/5000, dispersion: 2.55 A per pixel) and red grat- 
ing (600/7500, dispersion: 1.28 A per pixel) (see Lagattuta et al. 2010). 



the two measurements. For most lens galaxies, the agreement is 
good within the 68% confidence level error-bars measured for 
the photometric redshifts. Thus, these new photometric redshifts 
for the lenses are reliable estimates of their spectroscopic red- 
shifts. We notice that three systems show a strongly deviant pho- 
tometric redshift (COSMOS 0254+1430: z^^O.46+0.01 and 
z^ c =0.417+0.001; J095930.93+023427.7: z phot =lA0+™ 2 and 
z^ c =0.892+0.001; COSMOS 0018+3845: z^=0.44+0.01 
and z spec =0. 9755+0. 0003). In those three cases, the photome- 
try of the lens galaxy obtained from ground-based data is con- 
taminated by the close and bright images of the source, and the 
photometric redshift determination is consequently biased. 

2.3. Stellar masses of the lens galaxies 

We have derived a first set of stellar m asses using the Bayesian 
code described in iBundy et"aTI d2006). In brief, a data couple, 
made of the observed galaxy spectral energy distribution (SED) 
and its r edshift, is referenced to a grid of models constructed 
using the lBruzual & Charlol (120031) synthesis code. The grid in- 
cludes models that vary in age, star formation history, dust con- 
tent and metallicity. At each grid point, the probability that the 
observed SED fits the model is calculated: the corresponding 
stellar mass and the stellar mass to K-band luminosity ratio are 
stored. By minimizing over all parameters in the grid, the stel- 
lar mass probability distribution is obtained. The median of this 
distribution is taken as the stellar mass estimate, while the width 
corresponding to 68% of the probability distribution encodes the 
error bar resulting from degeneracies and uncertainties in the 
model parameter space. The final error bar on the stellar mass 
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Table 2. Position and redshift of galaxies in the FORS1 
field around COSMOS 0012+2015, COSMOS 0013+2249, and 
COSMOS 0049+5128. The error on the redshifts is +0.001. 



RA 


DEC 


z 


COSMOS 0012+2015 


150.01554 


2.3105342 


0.492 


150.05263 


2.3377222 


0.378 


150.07391 


2.3256596 


0.222 


150.07463 


2.3485102 


0.217 


150.07239 


2.3679651 


0.344 


150.08564 


2.3641732 


0.340 


150.10210 


2.3524922 


0.220 


150.11477 


2.3328146 


0.373 


COSMOS 0013+2249 


150.10193 


2.3900368 


0.218 


150.09001 


2.3836205 


0.530 


150.08210 


2.3904710 


0.351 


150.07843 


2.3895686 


0.352 


150.07239 


2.3679506 


0.347 


150.06510 


2.3855444 


0.343 


150.05808 


2.3804333 


0.346 


150.05235 


2.3837337 


0.350 


150.03149 


2.3745112 


0.220 


150.02774 


2.3738948 


0.222 


150.00332 


2.4075260 


0.347 


COSMOS 0049+5128 


150.23256 


1.8862956 


0.673 


150.23365 


1.8704970 


0.283 


150.21829 


1.8847159 


0.621 


150.20779 


1.8823028 


0.026 


150.20824 


1.8753444 


1.146 


150.20526 


1.8578028 


0.337 


150.19758 


1.8606685 


0.267 


150.18165 


1.8552841 


0.168 



also includes the K-band photometric uncertainty as well as the 
expected error on the luminosity distance that results from the 
uncertainty on the photometric redshift. The stellar masses and 
uncertainties are provided in Tableland plotted in Fig. [4j 

From the ACS-F814W images (Koekemoer et al. 2007) we 
infer the contamination due to the light coming from the back- 
ground lensed images in the stellar mass measurements. In an 
aperture of 3" radius, the lensed images are 3 to 7 magnitudes 
fainter than their respective lensing galaxy: in all cases this con- 
tamination generates an error on the stellar mass < 0.01 dex. 

In order to appreciate the reliability of our estimates of the 
stellar mass (and uncertainties), we have retrieved among the 
3.6yu IRAC catalog of the Spitzer-COSMOS survey (Sanders et 
al. 2007) 18 of the 20 strong lenses; the other two lens galaxies 
were not detected with Spitzer. For these lenses we have com- 
pared our estimates of the stellar masses to those inferred by 
Ilbert et al. (2O1C0): they are in reasonable agreement (see Fig. 
|4]). The typical difference between the two mass estimators (-0.2 
dex) is widely discussed in Ilbert et al. (2010). Both distributions 
show a tendency of increasing lens stellar mass with redshift. 
This feature will be discussed in more detail in the discussion 

(§EB. 



3. The environment of the lens galaxies 

Having now more accurate redshift measurements for the lens 
galaxies and the galaxies in the COSMOS field, we can under- 



2 The COSMOS cat alogs 
http ://irsa.ipac .caltech.edu/data/cosmos 
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Fig. 3. Photometric redshift, z p h t, versus spectroscopic redshift, 
z spec , for 16 of the lens galaxy candidates. The solid line features 
Zphot=Zspec- Errors bars on the value of z p h t relate to the 68% 
confidence level. 
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Fig. 4. The stellar masses (in M©) of the lens galaxies as a func- 
tion of redshift (spectroscopic when available, photometric oth- 
erwise). Black crosses: stellar masses inferred using the method 
described in Bundy et al. 2006. Grey stars: stellar masses in- 
ferred for the IRAC sources in Ilbert et al. 2010. Black and 
grey lines: the respective least square fits, with a slope of 
+0.8±0.3(black line) and +0.7+0.3 (grey line). 



take a study of the local environments of lens galaxies. Indeed, 
to understand whether high redshift lensing galaxies (z>0.3) are 
representative of the ETG population at their redshift, as the low 
redshift lens galaxies are (T09, Auger 2008), we study their en- 
vironment in comparison to the environment of a population of 
non-lensing ETGs. Our first analysis of the lens environment 
with regard to LSS (Paper II) led to the conclusion that lens 
galaxies are indeed evolving in the same environment than their 
parent population. Lets now zoom into a more local environment 
for the lens, using neighbor density measurements as defined and 
used by T09 for the SLACS sample. 
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Table 3. Summary of the properties of the lens galaxy environments (top to bottom: from low to high redshift). 



Name 


z 


i + + 8i + 


<2io>f 


Rio 


^1 
<Di> t 


Ni 


M* 






lens/source 




kpc 






lo g 10()M Q 




COSMOS 0049+5128 


0.337+0.001 


20.32+0.01 


1.3+1.9 


1013 


1.0+0.9 


10 


10.72+0.08 


63 


COSMOS 5947+4752 


0.345+0.001 


20.36+0.01 


1.9+1.5 


947 


1.6+0.7 


13 


10.65+0.22 


59 


COSMOS 0013+2249 


0.346+0.001 


19.63+0.01 


6.2+0.2 


602 


2.3+0.5 


16 


11.29+0.10 


23 


COSMOS 0056+1226 


0.361+0.001/0.808+0.001 


19.77+0.01 


1.4+1.2 


996 


1.1+0.7 


11 


11.26+0.13 


42 


COSMOS 5857+5949 


0.372+0.001 


20.55+0.01 


0.9+4.5 


1348 


0.4+2.2 


5 


10.81+0.09 


69 


COSMOS 0012+2015 


0.378+0.001 


19.95+0.01 


1.9+1.1 


833 


1.1+0.7 


12 


11.14+0.09 


40 


COSMOS 0254+1430 


0.417+0.001/0.779+0.001 


19.76+0.01 


0.5+3.5 


1767 


0.6+1.1 


5 


11.53+0.17 


28 


COSMOS 5921+0638 


0.551+0.001/3.14+0.05 


20.82+0.01 


0.6+4.1 


1939 


0.4+1.6 


3 


10.79+0.10 


51 


COSMOS 0216+2955 


0.608+0.001 


20.71+0.01 


4.6+0.3 


665 


2.7+0.3 


19 


11.38+0.07 


29 


COSMOS 0038+4133 


0.738+0.001 


21.03+0.01 


5.5+0.9 


585 


1.9+0.5 


17 


11.62+0.03 


30 


COSMOS 0124+5121 


0.856+0.001 


22.61+0.02 


2.9+1.0 


777 


2.1+0.4 


15 


10.78+0.08 


140 


COSMOS 0047+5023 


0.87+0.01 


21.35+0.01 


2.3+4.8 


851 


1.3+0.7 


11 


11.70+0.04 


22 


J100140. 12+020040.9 


0.879+0.001 


22.01+0.01 


0.6+4.7 


1931 


0.3+2.6 


2 


11.16+0.04 


94 


COSMOS 5941+3628 


0.88+0.01 


21.60+0.01 


2.5+0.4 


946 


2.0+0.3 


12 


11.54+0.06 


39 


J095930.93+023427.7 


0.892+0.001 


21.98+0.01 


1.2+2.4 


1131 


0.8+0.9 


7 


11.02+0.05 


92 


COSMOS 0211 + 1139 


0.92+0.01 


21.65+0.01 


1.8+0.6 


948 


1.6+0.4 


13 


11.76+0.04 


31 


COSMOS 0050+4901 


0.960+0.001 


22.01+0.01 


0.8+4.9 


1735 


0.6+1.5 


4 


11.64+0.03 


43 


COSMOS 0227+0451 


0.975+0.001 


22.09+0.01 


0.5+6.8 


2149 


1.0+0.9 


7 


11.39+0.04 


58 


COSMOS 0018+3845 


0.9755+0.0003/3.96+0.02 


23.43+0.03 


1.0+1.3 


1799 


0.3+1.2 


5 


10.73+0.06 


119 


COSMOS 5914+1219 


1.13+0.10 


23.00+0.02 


0.9+0.3 


1832 


0.7+1.2 


3 


11.39+0.05 


45 



Notes. (1) Lens name. (2) Redshifts of the lens ing g alaxies, followed by that of the source when known. (3) Subaru i + magnitude. (4) and (6) 
Number density ratios measured as explained in § 13.11 (5) Radius encompassing the 9 th closest neighbor. (7) Total number of galaxies within IMpc 
from the lens. (8) Lens galaxy stellar mass. (9) Number of twins in the photometric catalog of Ilbert et al. (2010). 



3.1. Projected number of neighbors 
3.1 .1 . The neighbor number densities 

In T09, the first estimator is the projected number density of 
galaxies, £io, inside a circle with radius Rio, equal to the pro- 
jected distance to the ninth closest neighbor of the lens. The 
neighbor galaxies are defined such as their magnitude is : i + < 
z'^+3, and the upper or lower bound, z, of their photometric red- 
shift is : Zlens ~ Aliens < Z < Ziens + dziens where Szi ens is a function 
of the lens magnitude and redshift, and lies between 0.012 and 
0.020 for our sample (see Ilbert et al. 2009). For the two lenses 
with incorrect photometric redshifts, we have used their spectro- 
scopic redshift instead and for the search of neighbors, the Szuns 
values associated with their spectroscopic redshift in the photo- 
metric catalog; this ensures that we are using the same search 
box for lens with spectroscopic and photometric redshifts. 

The search for neighbors is performed in the photometric 
redshift catalog of Ilbert et al. (2009) which contains about 1 
million objects. The limiting magnitude chosen for the neigh- 
bors ensures that the scaling in brightness (thus in mass) be- 
tween a lens galaxy and its neighbors is the same for all lenses, 
independently of its magnitude. For the faintest lens in our sam- 
ple, this upper limit corresponds to the completeness limit of the 
COSMOS photometric catalog (depth: i + ~ 26.2, Capak et al. 
2007a). The range chosen for the redshift ensures that, given the 
error bars on the neighbor photometric redshifts, they are indeed 
at the same redshift as the lens galaxies (in the 68% confidence 
limit error bars). 

As in T09, we also computed the projected number of galax- 
ies, D\, inside a circle with radius 1 Mpc, at the lens redshift. To 
select the neighbors, we have used the same limiting magnitudes 
and redshift ranges as for the first estimator. 

In order to get comparative measurements for a population 
of galaxies which are not identified as lenses but have the same 
morphology (early photometric type), similar brightness and 
redshift properties, we define a control sample of "twins" for 



each lens galaxy. Those twins have the upper or lower bounds, 
^twini of their photometric redshift such as: ziens - Sz < Ztwin < 
Ziens + <fe with (fe=0.05, and the upper or lower bounds, l tW i n , of 
their magnitude such as : \\ ens - SI < l tW i n < hens ~ SI, with 
£1=0.05 mag. The dispersion authorized around the value of the 
lens redshift and lens magnitude ensures that the sample of twins 
is large enough for the results to be statistically reliable (between 
22 and 140 twins, see Table [3]). In addition, the range in redshift 
is small enough to avoid possible biases due to galaxy evolu- 
tion and the range in magnitude ensures that the flux difference 
between a lens and its twins remains lower than 5%. We then 
measure the average number density of neighbors < Zio > t and 
the average projected number of galaxies in a circle with radius 
1 Mpc: < D\ > t for the control "twin" sample. 

We notice that, for an obvious lack of spectroscopic infor- 
mation, the twins cannot be selected according to their velocity 
dispersion as it is the case in the analysis by T09. Yet our match 
in magnitude and redshift mimics, up to a certain level of accu- 
racy, a match in mass between the lens galaxies and their twins. 

3.1 .2. Error measurements on the neighbor number 
densities 

In building the photometric catalog of galaxies, masks were nec- 
essary to hide bright objects disturbing the galaxy extraction 
procedure (see Capak et al. 2007a). From these masks we have 
measured, at each lens and twin location, the surface that was 
hidden for the extraction. The hidden surface is typically of a 
few percents of the total surface used to determine Eio or D\. It 
only gives a positive contribution to the total error budget on the 
neighbor density estimates: indeed, the surface covered by the 
circle of radius Rio would be smaller if the ninth closest galaxy 
was behind the mask, and if one or more galaxies are hidden be- 
hind the mask, the number of galaxies in a circle of radius D\ 
would be larger. Therefore, for a fixed number of galaxies, the 
surface has to be smaller. In summary, for a lens galaxy, the to- 
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Fig. 5. The distribution of 



(top plot) and 



(bot- 



<s 10 > f <D 1>f 
torn plot) as a function of the lens redshift, for the COSMOS 
lenses (red stars) and the SLACS lenses (black crosses). In each 
plot, the lower panel displays a re-binned version of the up- 
per plot using the the median of the ratios in different bins 
(SLACS bins: z<0.2 (32 galaxies), 0.2<z<0.3 (23 galaxies), 
0.3<z<0.5 (9 galaxies); COSMOS bins: 0.30<z<0.42 (7 galax- 
ies), 0.42<z<0.89 (7 galaxies), 0.89<z<1.30 (6 galaxies). 



tal error bar on the neighbor galaxy density is the quadratic sum 
of the Poisson error and the error due to the use of masks on 
the surface encompassing the neighbors. For its average twin, 
the error bar is the quadratic sum of the Poisson error, the error 
due to the use of masks and the dispersion for the population of 
twins. The Poisson error largely dominates the error budget for 



the lens neighbor number densities, while the dispersion largely 
dominates the error budget for the twins. 

3.1.3. Results 

The values of Eio, D\, < Eio > t and < D\ > t depend on the 
redshift bin and magnitude limit. To circumvent these effects, 



we only interpret the results in term of ratios: <£l 10 >f and <D l . 
These ratios quantify the richness of the environment of the lens 
galaxies in comparison to the environment of their respective 
twins. They are provided in Table [3] The average distance en- 
compassing the 9 closest galaxies to the lens is <Rio >=1.2 Mpc 
with a l<x standard deviation of 0.5 Mpc, hence comparable to 
the radius of 1 Mpc used to calculate D\. In Fig. we have dis- 



played the ratios 



and for the COSMOS lenses and 



<S 10 > f " AAV * <£>i> f 

for the SLACS lenses (using data from T09), as a function of 
redshift. A-priori, the two distributions cannot be directly com- 
pared as they may have different normalization factors (due to 
their different lens selection function, source redshift distribution 
and survey sensitivity). Yet, we notice that in the redshift range 
common to the SLACS and COSMOS samples (z~[0.33,0.50]), 



the levels of the distributions for 



and 



are similar, 



<Sio> f <Di> t 

suggesting that the normalizations are not drastically different. 
However, in order to avoid any misinterpretation due to intrinsic 
sample differences, the SLACS sample and the COSMOS sam- 
ple distributions are analyzed separately. 

The number density ratios displayed in Fig. [5] indicate 
whether the lens galaxies reside in a typical environment in 
comparison to the twin galaxy population. Would this be true, 
the number density ratios should be of the order of 1 at every 
redshift. This is indeed the case in average for the lenses in the 
redshift range studied here. 



3.2. X-ray and optically detected galaxy groups and clusters 

The validity of the estimators used to characterize the environ- 
ments of lens galaxies in § 13.1 l ean be checked using the distribu- 
tion of groups and clusters in the field based on X-ray and spec- 
troscopic data analysis. Indeed, we can access two supplemen- 
tary pieces of information: (a) the distribution of X-ray emitting 
gas in the COSMOS field (XMM-Newton survey: Hasinger et al. 
2007 and C-COSMOS with Chandra: Elvis et al. 2009), which 
traces galaxy groups and clusters in the field (Finoguenov et al. 
2007 and in prep.), and (b) the identification of optical groups or 
clusters (Knobel et al. 2009) from the z-COSMOS spectroscopic 
program (Lilly et al. 2007). 

Moreover, as the presence of groups in the direction to the 
source generates a different lensing potential in comparison to 
that of a single lens galaxy, it is important, in preparation of the 
mass modeling in § [4] to know: 1) If the lens galaxies are group 
members. If this is indeed the case, the group has to be modeled 
as an individual potential in the lens modeling. 2) What contri- 
bution in terms of shear is produced by the groups and clusters 
towards and around the lens galaxy. Most of the time this con- 
tribution is simply represented by an "external shear" whose pa- 
rameters are optimized while adjusting the lens model. In this 
paper, we rather intent to measure independently the external 
shear using the rich dataset that covers the COSMOS field. 

In the framework of our lensing analysis, each of the group 
catalogs (X-ray selected versus spectroscopically selected) has 
certain advantages and disadvantages. The entire COSMOS re- 
gion has been mapped through 54 overlapping XMM-Newton 
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Table 4. Summary of the contribution of the environment to the lens galaxies. 



Name 


P 


N 


K 




7 




PA r 

deg 


Cluster 

ID (z±5z, loglO(M 20 o), c, r 200 ) 


Dist. 

kpc 


COSMOS 0049+5128 


0.33 


13 


0.006+0.001^ 


TK50 — 


0.008+0.003+° 


T5UI — 


-38+12+° 


- 


- 


COSMOS 5947+4752 


0.04 


3 


0.016+0.003™ 


88A 


0.021+0.005™ 


884 


0+2+° 


- 


- 


COSMOS 0013+2249 


0.33 


26 


0.023+0.003™ 


M 


0.025+0.005™ 


88? 


8+10+} 


#173 (0.348+0.005,13.44,4.18,0.55) 


0.5 


COSMOS 0056+1226 


0.76 


29 


0.011+0.003™ 


888 


0.003+0.004™ 


888 


-18+30+° 


#133 (0.360+0.005,13.28,4.32,0.48) 


0.3 


COSMOS 5857+5949 


0.58 


18 


0.019+0.005™ 


W 


0.033+0.010™ 


88§ 


-41+8™ 


- 


- 


COSMOS 0012+2015 


0.39 


31 


0.019+0.002™ 


ill!.] 


0.009+0.006™ 


M 


-24+19+!! 


- 


- 


COSMOS 0254+1430 





2 


0.002+0.000™ 


883 


0.006+0.000™ 


ffi 


-31+3™ 


- 


- 


COSMOS 5921+0638 


0.61 


24 


0.019+0.004™ 


888 


0.009+0.012™ 


888 


-37+31™ 


- 


- 


COSMOS 0216+2955 


0.06 


5 


0.025+0.003^j 


88S 
885 


0.024+0.004+° 


888 


36+5+^ 


#221 (0.600+0.005,13.56,3.82,0.55) 


0.2 


COSMOS 0038+4133 


0.12 


4 


0.004+0.000™ 


0.014+0.002™ 


884 


" 4±4 -o 


- 


- 


COSMOS 0124+5121 


0.20 


2 


0.002+0.000™ 


881 
ijjH 

m 


0.007+0.002™ 


888 


9±3 -R 


- 


- 


COSMOS 0047+5023 


0.06 


4 


0.012+0.002™ 


0.016+0.003™ 


883 
88? 


-26+6™ 


- 


- 


n AAi /in 1 1 i ninr\Ar\ o 


U.ZU 




U.UUo+U.UUl q 




-28+7. 






COSMOS 5941+3628 


0.07 


7 


0.024+0.003+° 


w 

885 


0.007+0.004+° 


m 


22+15+? 






J095930.93+023427.7 


0.27 


4 


0.006+0.001 + ° 


0.016+0.004+° 


-45+1 


#288 (0.697+0.006,13.63,3.73,0.56) 


102 


COSMOS 0211 + 1139 





3 


0.004+0.000^ 


881 


0.022+0.005+° 


884 
887 
881 


43+3 + ? 






COSMOS 0050+4901 


0.25 


2 


0.001+0.000 + ° 


88! 


0.006+0.001 +° 


-16+4+5 






COSMOS 0227+0451 





4 


0.002+0.000+° 


88? 


0.005+0.001 +° 


-19+6^ 


# 101 (0.98+0.01,13.70,3.65,0.53) 


0.5 


COSMOS 0018+3845 


0.38 


15 


0.037+0.007+° 


88A 


0.082+0.016+° 


888 


13+4+0 






COSMOS 5914+1219 


0.25 


17 


0.001+0.000+° 


888 

000 


0.006+0.001+° 


888 

000 


-19+5i 







Notes. (1) Lens name. (2) Probability that galaxies in a 2'x2' field around the lens galaxy are included in the z-COSMOS catalog. (3) The number 
of groups involved in the calculation of the external shear (within 5' projected radius). (4, 5, 6) Convergence, shear strength and shear direction 
due to the groups and clusters around the lens. The first set of error comes from the uncertainty on the group mass. The second set of error are 
due to the uncertainty on the source redshift, chosen to be +0.5 for sources with unknown redshifts and the error on the source redshift estimation 
from Table [3] when it has been measured. (7) The ID of the group or cluster with an impact parameter smaller than its radius (r 2 oo) in Leauthaud et 
al. (2010). In parenthesis are its redshift, mass (in solar mass), concentration and radius (in Mpc). The error on the cluster mass is 40% of the total 
mass (Vikhlinin et al. 2009). (8) Projected distance in kpc between the lens and the group/cluster centre (at the lens redshift). 



pointings while additional Chandra observations have mapped 
the central region (0.9 square degrees). A composite XMM- 
Newton and Chandra mosaic has been used to detect and mea- 
sure the fluxes of groups and clusters to a 4<x detection limit of 
1.0 x 10" 15 erg cm" 2 s" 1 over 96% of the ACS field. The gen- 
eral data reduction process can be found in Finoguenov et al. 
(2007, 2009) and the cluster catalog and weak-lensing analysis 
are presented in Leauthaud et al. (2010). Therefore, the cata- 
log of X-ray selected groups and clusters is homogeneous and 
covers the entire COSMOS field. However, it is affected by the 
sensitivity threshold of the X-ray survey, which is different from 
the sensitivity cut of the optical dataset. Therefore, groups of 
modest mass could be missing from the X-ray selected catalog, 
especially around galaxies at high redshift. 

Conversely, the catalog of optically selected groups in 
COSMOS has the advantage of spanning the same brightness 
range as the optical imaging dataset from which the strong lens 
candidates and the arcs have been extracted. Even though the 
faintest galaxies seen on deep images will remain undetected 
in the shallower spectroscopic z-COSMOS survey, groups can 
still be traced and assessed from their central brightest galax- 
ies. Then, the parameters of the galaxy groups can be derived by 
comparing these detections to a catalog of mock-groups subject 
to the same detection criteria (see Knobel et al. 2009 for an ex- 
tensive description of the group catalog). The first data release 
of the z-COSMOS survey (Lilly et al. 2007) covers a fraction 
of the full COSMOS field, which leaves some of our COSMOS 
lens candidates outside the coverage of the group catalog. 

Let us establish now whether the lens galaxies are group or 
cluster members (§ 13.2.1b and measure the external shear con- 
tribution of the groups and clusters at the lens galaxy position 
(§[323]). 



3.2.1 . Lens galaxy group members 

We cross-correlate the X-ray cluster catalog (Leauthaud et al. 
2010) with the sample of strong lenses. We define as a group 
member, a lens galaxy located within r2oo (radius where the mat- 
ter density is 200 times the critical density) of the group cen- 
ter, and with redshift identical to the group redshift (within error 
bars). We have found 4 matches and there is a galaxy cluster de- 
tected at a lower-redshift in direction to J095930.93+023427.7 
(see Table 

Using the group catalog built from the optical dataset, we do 
not identify any new groups associated in redshift and observa- 
tional plane with lens galaxies than those already identified using 
the X-ray catalog. In Tabled we give the probability that galax- 
ies in a 2'x2' field around the lens galaxy are included in the 
z-COSMOS catalog. This figure establishes roughly the "com- 
pleteness" of the survey at every lens location and tells us if the 
z-COSMOS group catalog can be used to characterize the envi- 
ronment of the lens. 

3.2.2. Comparison: galaxy group members and projected 
number density of neighbors 

In principle, lens galaxies which are group members 
are expected to have large neighbor number densities. 
This should be visible in the ratio measured in § 13.11 
This is indeed the case with the first estimator for 
the lens in COSMOS 0013+2249 (-f^ =6.2+0.2) and 

COSMOS 0216+2955 (-f±V =4.6+0.3), and it is unclear for 
COSMOS 0056+1226 and COSMOS 0227+0451 because 
of the large error bar on the measured density ratios. With 
the second estimator, the ratio are in average lower than with 



7 



C. Faure et al.: Environment and mass properties of COSMOS lens galaxies. 



^fV> an d it i s therefore more difficult to identify correctly an 
over-dense field. On the contrary, COSMOS 0038+4133 shows 
a large density ratio =5.5+0.9) without being associated 

to any known galaxy group in X-ray or optic. 



3.2.3. The external shear due to clusters and groups 

We estimate the shear and convergence produced by the groups 
(either X-ray or optical) detected around the lens galaxies in or- 
der to characterize the contribution of the environment to the 
total lens potential in future lens models C§ SJ. To do so we fol- 
lowed the method described in Paper III. In short, we compute 
the convergence and shear produced by all groups closer in pro- 
jection than a given radius (5', as in Paper III), and at any red- 
shift up to that of the source. We assume that the mass profile 
of every group follows a truncated isothermal sphere (TIS). The 
choice of this profile is motivated by the fact that the shear and 
convergence produced by an isothermal sphere are easy to com- 
pute (Keeton 2003, Momchevaet al. 2006). Moreover, as we will 
consider only groups and clusters which do not cross the line-of- 
sight to the source, estimations of their total masses rather than 
of their mass distributions are sufficient for the accuracy of the 
result. We selected a truncated profile to avoid to give unrealisti- 
cally too much weight to the most distant groups. 

The three-dimensional density distribution of the TIS can be 
written: 



1 1 



p oc 



r 2 r 2 + r 2 



(1) 



where r c is the truncation radius. The convergence, k, and shear, 
y, produced by this profile are respectively: 



1 



1 



2 ^ Vr 2 V r ? + r 2 



V^ 2 ^ 



2r c { Vr 2 + r 2 _ ^ 



ri + r L 



(2) 



(3) 



where b is the impact parameter of the TIS, which relates to the 
Singular Isothermal Sphere (SIS) impact parameter bsis as: 



bsis 
b 



1 + 



bsis 



1 + 



bsis 



(4) 



oo these quantities match those of the 



In the limit where r c 
SIS. 

To compute the external shear produced by the groups and 
clusters we proceed as follows. Whenever P<0.3 in Table|4] (col- 
umn 2), we only consider the catalog of X-ray detected groups 
and clusters. The groups mass potentials are modeled as TIS 
using the value M200 and r2oo (used as the truncature radius). 
The error on the X-ray mass comes from the scatter in the re- 
lation used to derive the mass from the luminosity (Leautaud et 
al. 2010): it is of the order of +20% of M 20 o (Vikhlinin et al. 
2009). Whenever P>0.3 in Tabled (column 2), we also consider 
the optical group catalog in addition to the X-ray group catalog. 
We first correlate the optical and X-ray group catalogs to iden- 
tify and remove optical groups that might already be accounted 
as an X-ray group. Then, we model as TIS the remaining opti- 
cal groups using the "mock" virial mass, M VI >, and the "mock" 
virial radius, r vir (used as the truncature radius). The upper and 
lower error on the virial mass are of M VI > (Knobel et al. 



2009). These "mock" quantities are the theoretical values asso- 
ciated with the detected groups when subjecting a mock sample 
of groups to the same selection function and same survey criteria 
as the observations (see Knobel et al. 2009 for details). 

The external shear and associated convergence are calculated 
for each group individually. They are then re- scaled to the red- 
shift of the lens using the scaling relation given in Momcheva 
et al. (2006). For a given lens, external shear and convergence 
are summed up following the procedure described in Keeton 
(2003) and Momcheva et al. (2006). The results are summarized 
in Table ffl A first set of errors on the shear and on the conver- 
gence results from the propagation of the group mass errors. 

As mentioned already, in the special case where a galaxy 
group or cluster has an impact parameter smaller than r2oo or r vir 
(in 5 cases, see Table |4]), the shear and convergence calculated 
under this simple approximation are incorrect (Keeton 2003). 
Hence, we have systematically removed these groups when com- 
puting the external shear. Instead, they will have to be accounted 
for additional lens potential when performing the lens modeling. 

Regarding the source redshift, it is either known spectro- 
scopically for some lenses (see Table 0), or assumed to be at 
z s =2 for lenses with z/<l, or at z^=3 for lenses with z/>l. The 
convergence and shear contributions depend on the number of 
groups taken into account, which in turn may depend on the 
source redshift and on the cut-radius considered. Hence, an er- 
ror on the source redshift generates an error on the external 
shear parameter. We have estimated this error assuming an un- 
certainty fe=+0.5 on the source redshift when unknown, and us- 
ing the error on the source redshift in Table [3] when measured 
(see Table |4]). To estimate the error introduced by an arbitrary 
cut-radius at 5', we have also probed two other radii (7' and 
10'). We find that the incidence of the radius selection on the fi- 
nal lensing contribution is 6k ~ Sy -0.001 and is negligible on 
the orientation of the shear ( 5¥A y <0.5°). Then we calculate the 
total error on the external shear as the sum in quadrature of the 
error produced by the group mass uncertainty, the source redshift 
uncertainty and the error produced by the choice of aperture. 

We have also analyzed the error on the external shear pa- 
rameters generated when using a single catalog of groups (X- 
ray) instead of the combination of the two catalogs (X-ray and 
optical). The comparison is possible for 7 systems. We notice 
that, while the convergence is different when using a single cat- 
alog instead of two, the external shear strength and direction are 
in agreement within their respective error bars. This means that 
our calculation realistically associates the largest source of shear 
with the more massive groups and clusters. It also means that the 
probable incompleteness of our catalogs in low mass groups has 
a minimal impact on the shear measurement: what matters is the 
completeness of the catalog in large mass groups and clusters. 
However, the technique is limited for high redshift sources, as 
there is no X-ray cluster detected above z=1.3 in the COSMOS 
field. 



4. The strong lens modeling 

In this section, we focus on the 12 lens galaxies offering 
the largest number of observational constraints: the triple and 
quadruple image systems and the Einstein rings. 



4.1. Lens galaxy light profiles 
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Table 5. Lens galaxy luminous profile. 



Name 



COSMOS 0049+5128 
COSMOS 5947+4752 
COSMOS 5921+0638 
COSMOS 0038+4133 
COSMOS 0124+5121 
COSMOS 0047+5023 
J100140. 12+020040.9 
COSMOS 5941+3628 
J095930.93+023427.7 
COSMOS 0050+4901 
COSMOS 0018+3845 
COSMOS 5914+1219 



H 99+0.00 



Sersic parameters 
PA R e R e 

deg kpc " 



n 



1.19+0-01 
1.38+ 




Notes. (1) Lens name. Parameters of the Sersic light profile fit: (2) ellipticity, (3) position angle, (4) effective radius in kpc and (5) in arc-seconds 
and (6) index, with associated error bars (from GIM2D, 68% confidence limit) . 



We have re-computed the 2-dimensional fit of the galaxy 
surface brightness distribution of the COSMOS lenses using 
GIM2D (Simard 1998, Marleau & Simard 1998) in order to in- 
clude error bars that were not presented in Paper I. For that pur- 
pose, we adopt, as in Paper I, a Sersic bulge plus an exponen- 
tial disk parametrization to describe the two-dimensional sur- 
face brightness distribution of the lens galaxy light profile. The 
Sersic profile is parametrized by means of the total flux in the 
bulge, the Sersic index, n, the bulge ellipticity, 6=l-b/a, the posi- 
tion angle of the bulge, PA, and the effective radius of the bulge, 
R e . The exponential profile depends on the photometric disc to- 
tal flux, the disc scale-length, the disc position angle and the disc 
inclination. The software gives the best fitting values for all of 
these parameters. The parameters of the Sersic bulges are sum- 
marized in Table [5l The error bars correspond to the 68% con- 
fidence level. For most systems, the results are consistent with 
the surface brightness parameters measured in Paper I. But for 
others, such as COSMOS 5921+0638, the best fit parameters are 
different in Paper I, Paper III and here. Indeed, the presence of a 
ring or bright arcs close to the lens galaxy center makes difficult 
to produce a robust fit of the lens galaxy surface brightness den- 
sity profile; this remains the case even when more complex fit- 
ting and deconvolution methods are used (see Chantry & Magain 
2007). 



The relative image positions to the lens galaxies are the main 
constraints for the lens models. For images which are point like 
objects, the determination of their positions only depends on the 
image resolution. This is the case for COSMOS 5921+0638. 
For this system the error on the relative position is 0.014" (see 
Paper III). For the other multiple images systems, we determine 
the position of the the brightest peak in each image. For the- 
ses systems, the error on the relative position of the images is 
typically 0.05". For the perfect rings that do not display any 
peak we place the image arbitrarily around the ring, in a sym- 
metric way around the lens center and we assume that the er- 
ror on the relative position of the images is 0.05". In Table 
we provide the lens galaxy central coordinates as well as the 
position of the multiple images relative to the lens galaxy lo- 
cation used in the lens modeling. For J100140. 12+020040.9 
and J095930.93+023427.7, we have retrieved the image position 
from Jackson 2009. 



4.2. The mass models 

The purpose of the new lens modeling is: i) to measure the total 
mass of the galaxy within R#: M(<R#) and ii) to check if the 
environment contribution measured in § !3.2.3l is a good estimator 
of the external shear . 

To do so, we have used the Lenstool code (Kneib et al. 1993, 
Julio et al. 2007) to model the lens galaxy mass distributions. 
Lenstool allows a x 2 minimization of parametric mass models, 
either in the source or in the image plane. For higher accuracy, 
we use the image plane minimization algorithm. 

4.2.1. SIE+y 

We have performed a set of lens model in which the mass dis- 
tribution the lens galaxies follows that of a Singular Isothermal 
Ellipsoid (SIE). A SIE is defined by its position, velocity dis- 
persion, orientation and ellipticity. Assuming that mass follows 
light, we fix the SIE central position to that of the galaxy light 
profile (Table©. We chose its orientation to be that of the Sersic 
bulge light profile (Table [5]) with an additional +10° uncertainty 
to take into account for a possible misalignment between the 
luminous bulge and the DM halo (e.g. Kochaneck 2002). The 
higher boundary for the SIE ellipticity is set equal to that of the 
Sersic profile (Table 0; the lowest boundary is set to in order 
to take into account for a possibly shallower distribution of the 
DM halo in comparison to the luminous core observed in ETGs 
(e.g. Gavazzi et al. 2007). 

In addition, we model the contribution of the environment 
by an "external shear" which is parametrized by a shear strength 
and orientation. The group and cluster contrib utions to the shear 
at the lens galaxy locations are measured in § 13.2.31 We use the 
values of y and PA r from Table ffl and the corresponding error 
bars as priors, to build the lens models. 

Therefore, the total number of free parameters for the models 
is 5: 3 for the SIE (orientation, ellipticity, velocity dispersion) 
and 2 for the external shear (strength and direction). The number 
of observational constraints is 6 for the quads and rings and 4 
for the triple image lenses. The x 2 of me best lens models are 
reported in Table |7] (column 3). The best fit models (with^ 2 in 
column 3) are displayed in Figs. [7] to [TT] Whenever^ 2 » 1, 
we have performed a second lens model, changing the priors 
on the external shear to the following ones: y- [0.0,0.9], PA r =[- 
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Table 7. Best fit parameters for the lens models: SIE + shear. 



Name N 


x 1 


Ct 2 ',y,PA r ) 


cr v + (5cr v R/. 


R £ 


M(<R £ ) 


fDM(<R^) 








km s~ l " 

ai a +7 01 t+U.UU 


kpc 


70 n M o 

-7 a 1+0.11 


r» n /I +0.06 



COSMOS 5947+4752 R 1.6 

COSMOS 5921+0638 4 11.8 

COSMOS 0038+4133 4 6.3 

COSMOS 0124+5121 R 0.9 

COSMOS 0047+5023 4 1283.1 

J100140. 12+020040.9 4 7.3 

COSMOS 5941+3628 R 0.5 

J095930.93+023427.7 4 6.8 

COSMOS 0050+4901 4 545.4 

COSMOS 0018+3845 3 45.3 

COSMOS 5914+1219 3 215.9 



(2.0, 0.025, +20.5) 
(0.5,0.128, -76.3°) 

(10.4, 0.126,-15.10°) 
(2.3,0.053, +0.16°) 



(26.3, 0.097, +79.9°) 
(0.2, 0.103, +39.9°) 
(0.1,0.086, +71.2°) 



326+J 1 
189 + ° 
207^3 

383 -38 

959+43 

255 +s * 
386 ; § 
289! 2 
358!% 



2.33+] 
0.62+1 

o.6o + ; 

0.9l + _\ 
1.80+! 

o.8r; 

1.23+1 
0.79+1 

1.65+; 

1.32+1 
1.60 ' 



0-01 



11.1+ 
4.4! 
4.3+^ 
6.8+| 

13.5!? 

6.0+[ 
12.7! 
10.2!J 
12.8 



8.61 

1.10: 

3 55+^ 

-0.5: 

14.52+^! 
3.00+ i:ii 



6.80" 
2.89; 

11.88 
6.22" 

12.00" 



0-87 



0.95!| 
0.58!j 
-0.37!| 
0.84+; 
0.68+; 

0.6O!; 

0.62!; 

0.71!; 
0.92!; 

0.93!; 
0.79!; 



3 

0-04 



Notes. (1) Lens name. (2) Number of images or "R" if it is a complete Einstein ring. (3)x 2 of the best lens modeled obtained when using the priors 
on the external shear measured in § 13.2.31 (4) Favorite shear parameters (when the one measured in § !3.2.3l lead to^ 2 »1) and corresponding^- 2 . 
(5) Velocity dispersion of the SIE corresponding to the best model in (4) or in (3) if no other. (6) Corresponding Einstein radius in arc-second and 
(7) in kpc. (8) Mass of the lens galaxy in the Einstein radi us. (9) Fraction of DM in the Einstein radius. In column (5), (6), (7) and (8), the errors 
are due to the source redshift uncertainty (as explained in § 13.2.31 ). 



90°, 90°], hence letting the external shear parameters free. The^ 2 
of these second lens models are reported in Table [7] (column 4) 
and referred as x 2 ' • 

For J095930.93+023427.7, we have modeled the group in 
direction to the lens by a SIS which position is fixed to the posi- 
tion of the group in the X-ray catalog. Hence, the only parameter 
allowed to vary is the velocity dispersion of the profile. The best 
fit model in Table [7] is obtained for a group with velocity disper- 
sion: al roup =408+29 km s" 1 , if it was at the lens galaxy redshift. 
From Fig. [TT] we see that the fit is not perfect: with additional 
observational constraints such as the velocity dispersions of in- 
dividual group members and of the lens galaxy closest neighbor 
(to the East), one could perform a more detail model of the lens 
potential which would most probably improve the fit. 

We notice that systems with perfect Einstein rings are mostly 
satisfactorily modeled by a SIE plus the measured external 
shear (COSMOS 5947+4752 in Fig. COSMOS 0124+5121 
in Fig. HI COSMOS 5941+3628 in Fig. The reason is that, 
for those lenses, we have arbitrarily chosen the image positions: 
they are symmetrically distributed around the lens. In addition, 
the SIE ellipticity is allowed to be null. Hence we are artificially 
correctly fitting the image positions, whatever the external shear 
values are (as long as the shear strengths are not too large). For 
those systems, only the Einstein ring and associated mass are re- 
liable measurements in Tabled 

For COSMOS 0049+5128, the fit is not good 0t 2 =35.3). If we 
let the external shear free, the best fit shear parameters are dif- 
ferent than the one measured in § 13.2.31 If we subtract the shear 
vector given in Tabled to the best fit shear vector in Tabl^T] (col- 
umn 4), we find the direction pointing toward the galaxy closest 
to the lens (Galaxy 2 in Fig. [7]). Galaxy 2 is at a projected dis- 
tance ~13" to the lens galaxy. Would it be at the lens redshift, 
it would need to have a velocity dispersion (T ^ /ax >' 2 -i74 k m s -i 
to create the shear necessary to obtain;^ 2 ' while fixing the shear 
parameters to the measured values. 

It is the same for the other lenses: if fixed in the lens model, 
the external shear due to the groups leads to^ 2 » 1. If we set 
the shear parameters free, the best fit shear will point in a direc- 
tion which correspond to the vectorial summed orientation of the 
shear due to the groups and of the shear due to a secondary (and 
third in the case of COSMOS 0038+4133) galaxy. In Table [8] 



we have reported the projected distance between the lens galaxy 
and the secondary (third) galaxy, as well as the the velocity dis- 
persion of the second (third) galaxy that is needed to obtain x 1 ' 
when fixing the external shear parameters to the one due to the 
groups. In the case of COSMOS 0050+4901, the velocity disper- 
sion derived for the main lens and for the secondary galaxy are 
that of groups rather than that of galaxies: indeed, when we look 
at the image of the lens (Fig.[l0|), we see that the field is crowded 
with galaxies. However, the neighbor density ratio measured in 
§ [3] do not show any evidence for the presence of a structure at 
the lens redshift. Therefore if there is actually a group interven- 
ing in this system, it should be at a different redshift than that of 
the main lens galaxy. For the other systems, the velocity disper- 
sion associated with the secondary lens can be associated with a 
galaxy mass. For COSMOS 0047+5023, we see in Fig.[9]that the 
field around the lens is crowded with galaxies: an improved ver- 
sion of the present mass model should take them into account, 
preferentially using a measurement of their redshift and velocity 
dispersion. 

We conclude that our efforts at measuring the external shear 
were not in vain as, when using it in the lens model, we clearly 
identify the missing element : the best fit will point at the closest 
galaxy to the lens or indicate the realistic presence of a galaxy 
group in the line-of- sight. We come back to this result in § 15.41 

The Einstein radius of the lens galaxy and corresponding 
mass have been calculated for the best models and are displayed 
in Tabled 

4.2.2. The proportion of DM in the Einstein radius 

We integrate the galaxy light profile density up to the Einstein 
radius, and use the stellar mass of the lens galaxy, identifying 
the effective radius to the half stellar mass radius, to determine 
the stellar mass within the Einstein radius, M*(<R#). Doing so, 
we can compare it to the total mass within the Einstein radius, 
M(<R#), obtained during the lens modeling. This gives us a 
measurement of the lens galaxy projected DM fraction within 
the Einstein radius: /dm(<R^)= I' l^T) • These values are re- 
ported in Table |7J The errors come from the propagation of the 
uncertainties in M*, R e , M(<R#) and R#. 
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Table 6. Position of the lensing galaxy (in degree) and relative 
position (in arc-second) of the images used as constraints for the 
lens mass model. 



Lens 


RA 


DEC 


Image 


5RA 


6DEC 


COSMOS 0049+5128 


150.2052807 


1.8578028 


Al 


2.05 





A2 





2.05 


A3 





-2.05 


A4 


-2.05 





COSMOS 5947+4752 


149.94955 


2.7979502 


Al 


2.28 





A2 


-2.28 





A3 





2.28 


A4 





-2.28 


COSMOS 0038+4133 


150.1595 


2.6927351 


Al 


-0.43 


-0.33 


A2 


0.00 


-0.62 


A3 


0.04 


-0.61 


A4 


0.22 


0.71 


COSMOS 0124+5121 


150.3522 


1.8558994 


Al 





0.89 


A2 





-0.89 


A3 


0.89 





A4 


-0.89 





COSMOS 0047+5023 


150.19858 


1.8397919 


Al 


-1.42 


-1.15 


A2 


0.10 


-1.56 


A3 


2.08 


-0.36 


A4 


-0.05 


0.90 


COSMOS 5941+3628 


149.92209 


2.6080412 


Al 


0. 


1.24 


A2 


0. 


-1.24 


A3 


-1.24 





A4 


1.24 





COSMOS 0050+4901 


150.21104 


2.8171935 


Al 


1.15 


0.69 


A2 


-1.53 


1.10 


A3 


-1.53 


-0.80 


A4 


0.27 


1.88 


COSMOS 0018+3845 


150.07666 


2.6458333 


Al 


-0.39 


-1.15 


A2 


-1.22 


-0.33 


A3 


1.36 


0.25 


COSMOS 5914+1219 


149.81142 


2.2054236 


Al 


-0.93 


-0.45 


A2 


1.26 


1.43 


A3 


1.81 


0.55 



Table 8. Parameters for Galaxy 2. 



Name 


Distance 

" 


OalaxyT. 

km s~ l 


Fig. 


COSMOS 0049+5128 


12.8 


174 


HI 


COSMOS 5921+0638 


1.6 


70 


8 


COSMOS 0038+4133 


3.5 


195 


8 


// 


4.6 


195 


8 


COSMOS 0047+5023 


3.9 


264 


9 


J100140. 12+020040.9 


1.1 


94 


Q 


J095930.93+023427.7 


4.9 


275 


m 


COSMOS 0050+4901 


7.4 


366 


M 


COSMOS 0018+3845 


3.5 


113 


M 


COSMOS 5914+1219 


3.0 


209 


M 



Notes. (1) Lens Name. (2) Distance between the lens galaxy and 
Galaxy 2 (or Galaxy 3 in the case of COSMOS 0038+413 3). (3) 
Velocity dispersion of Galaxy 2 (or Galaxy 3) as explained in § 14.2.11 
(4) Figure number where the lens model is displayed and the Galaxy 2 
and 3 are labelled. 



For the 11 remaining systems, the fraction of DM varies 
between 0.58+0.00 (COSMOS 5921+0638) and 0.95^ 
(COSMOS 5947+4752). 

In Fig. [6] we have reported the evolution of foM (<R#) with 
the redshift. The evolution of f^M(<R^) with redshift is com- 
patible with constant or slightly decreasing. The interpretation 
of this tendency is linked to the evolution of the Einstein radius, 
Re versus effective radius, R e . Indeed, for a given galaxy if (1) 
Re ~ R e , the fraction of dark matter in the Einstein radius is 
expected to be "low" as, in this case, we are probing the region 
where the galaxy is baryon dominated: f^M(<R^) is then a 
lower limit of the total foM- On the contrary, if (2) R E » R e \ 
we are measuring a f^M(<R^) which is getting close to the total 
foM, and therefore it is expected to be larger, in average, than 
the one determined in case (1). This is indeed the case (see 
FigO top right panel): foM(<R^) increases slowly when Rs/Re 
increases. From the same Fig. [6] (bottom left panel) we see that 
the ratio Rn/Re increases quickly with redshift, meaning that 
as redshift grows we are measuring fDM(<R#) getting closer to 
the total foM of the galaxy. The fact that foM(<R^) is slightly 
decreasing or "at best" constant when the redshift grows suggest 
that the "total" fraction of dark matter is genuinely lower in high 
redshift lens galaxies, than in the low redshift lens galaxies. 



The negative DM fraction for COSMOS 0038+4133 may in- 
dicate that the stellar mass of the lens galaxy is overestimated. 
This is surprising as for this system the photometric and spectro- 
scopic redshifts agree, hence making us confident that the SED, 
hence the stellar mass, is reliable as well. Another possible ex- 
planation is that the light profile fit is not correct and that we 
measure too large an effective radius. This is indeed a possi- 
ble explanation as a bright ring surrounds the lens galaxy and 
may bias the determination of the light profile. Finally, it is also 
possible that the source redshift is lower than the one consid- 
ered here (zf in =1.5). Indeed, would z s =0.9, the total mass in the 
Einstein radius would be M(<R£)=4.88xl0 n M o , and the cor- 
responding DM fraction: f DM (<R^)=0.6. Because of all these 
interrogations, which requires further investigation, we do not 
keep this system for the rest of the analysis. 



Interestingly, we observe the same tendency if we consider 
the lens sample from Jiang and Kochanek (2007, JK07 here- 
after): they measure the DM fraction (assuming or not adiabatic 
compression in their galaxy models and Salpeter initial mass 
function (IMF)) in the Einstein radius of 22 galaxies spanning 
a lens redshift z/=[0.0808, 1.004]. If we look at the 18 galaxies 
with Re I R e <10 as in our sample (see Fig. [6]), we notice that, 
similarly to our sample: (i) f^M(<R^) decreases slightly with 
redshift with the same slope than for our sample (-0.18), (ii) the 
radius Rs/Re increases with redshift and (iii) f^M(<R^) slightly 
increases when Re/Rc increases. 

We add that, even if the two samples show similar behavior, we 
chose not to mix them for the display and slope calculations in 
order to avoid misinterpretations due to possible systematics af- 
fecting the measurement of the DM fraction in the two different 
methods (e.g. different IMFs). 
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Fig. 7. Mass models on top of ACS images for (top to bottom) COSMOS 0049+5128 and COSMOS 5947+4752. The left column 
displays the best mass model corresponding to the^ 2 given in Table [7] (column 4 or column 3 if no other). The right panel shows 
a zoom of the central region of the lens models, where the source can be seen. North is to the top and East to the left. Color code: 
the red circles are the observed images, which radii correspond to the position uncertainty used in the modeling. In orange are 
the images produced by the best lens model (in case of perfect fit, the orange and red crosses superimpose). The navy blue lines 
describe the potential. The caustic lines are in yellow an the critical lines are in cyan. The green ellipses show the position of the 
source as seen through the best mass model (one source for each image, when a good fit is reached the four sources are partially 
superimposed). We indicate the position of Galaxy 2 (see § 14.2. H and Table [8]). 



5. Discussion and conclusions 

From the analysis of the COSMOS strong lenses in the redshift 
range 0.34 to 1 . 13, we get three major results: (1) the lens galaxy 
stellar masses increases with the redshift, (2) the lens galaxy en- 
vironments are compatible with those of ETGs in this redshift 
range and (3) the DM fraction of the lens galaxies in the Einstein 
radius slightly decreases with the redshift, when the ratio be- 
tween the Einstein radius versus the effective radius strongly in- 
creases with redshift. Let us discuss each one of these trends. 



5.1. Increase of the lens stellar mass with redshift 

In § I2.31 we have measured that the lens galaxy stellar mass in- 
creases with redshift. Is this a result of the selection method used 
to built the lens sample? At z~0.4, one of our lenses has an effec- 



tive radius of -0.3" and at z~0.9, the smallest lens galaxy effec- 
tive radius is -0.2" (see Table [5]) . The effective radius is a lower 
limit for the Einstein radius. These "lower limit Einstein radii" 
correspond to lower limits in mass within the Einstein radius that 
are in our case: 10 10 9 M Q at z=0.4 and 10 110 M© at z=0.9, as- 
suming a source at z=1.5. Hence, for a given source redshift, the 
minimal lens galaxy mass should increase with the lens redshift 
in order to be visually detected; hence, most probably, the mini- 
mal stellar mass might increase as well with redshift. Following 
this reasoning, we expect the sample average stellar mass to in- 
crease with the redshift. While this reasoning is simplistic as the 
source redshifts are different for every lens, it still explains par- 
tially the increasing stellar mass with increasing redshift of the 
lens galaxies. So far, it is not clear whether this "detection bias" 
effect is the only reason for the lens galaxy stellar mass increase 
with redshift. This result should be investigated in dedicated nu- 
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Fig. 6. Evolution of the dark matter fraction and related parame- 
ters for two lens samples: COSMOS (crosses) and JK07 (stars). 
Top left panels: Evolution of the DM fraction in the Einstein 
radius with the redshift. The least deviation fits are displayed 
(solid lines), and follow the relations //)m(<R^)=-0.18xz+1.00 
for COSMOS and / DM (<R£)=-0.18xz+0.38 for JK07. Top right 
panels: Variation of the DM fraction as function of the ratio of 
the Einstein radius versus effective radius. The least deviation 
fits give the following relations: /dm(<R£)=+0.11x|£+0.41 

(COSMOS) and f DM (<R E )=+0.09x^+0.0& (JK07). Bottom 
panels: Evolution with redshift of the Einstein radius ver- 
sus effective radius. The least deviation fits give the follow- 
ing relations: §^=+3.4xz+0.8 (COSMOS) and §^=+4.3xz+0.7 
(JK07). For JK07, we have reported the results obtained for 
galaxy models without adiabatic compression. The slopes are 
similar for models with adiabatic compression. 



merical simulations (e.g. van de Ven et al. (2009), Mandelbaum 
et al. (2009)). 



5.2. Evolution of the lens galaxy environment 

Our study of the environments of COSMOS lens galaxies, via 
the projected galaxy number density in § I3.U indicates that the 
environment of lens galaxies is similar to that of ETGs across the 
whole range of redshifts tested here. This result extends the pre- 
vious measurements made in the redshift range [0.068; 0.513] 
with the SLACS sample (T09 and Auger 2008). It is particularly 
interesting as the COSMOS and SLACS samples have very dif- 
ferent selection criteria: hence this result is most probably a gen- 
uine characteristic of strong lens galaxies rather than a selection 
bias. 

We also notice that both neighbor density estimator tested 
here appear to be reliable estimators of the environment of galax- 



ies, as X-ray groups and clusters are actually detected around 
lenses with high < ^ >t and high tttV ratios. 



<Di> t 



5.3. The decreasing dark matter fraction with redshift 

We have measured the lens galaxy DM fraction in the Einstein 
radius by combining the total mass in the Einstein radius, the 
light density profile and the stellar mass in the galaxy. The pro- 
jected DM fraction decreases with redshift even though the ratio 
between the Einstein radius and the effective radius increases 
with the redshift (see Fig. [6]). A similar trend is seen in the lens 
galaxy dataset of JK07 which covers a similar redshift range. 
Using toy models and a ACDM cosmology and the SLACS 
lenses, Napolitano et al. (2010) and Tortora et al. (2010) first 
notice that, for a fixed stellar mass and age, the DM fraction is 
similar for high and low redshift galaxies. They also notice that 
the slope of the distribution of foM versus age is steeper than 
explained by their model and they invoke different scenarios to 
interpret the discrepancy (including adiabatic compression and 
different IMF for different galaxy ages). Another possible effect 
responsible for the DM fraction decrease with redshift is the in- 
creasing stellar density of ETGs with redshift as discussed by 
Bezanson et al. (2009). But to understand if the results are af- 
fected by one of this effect, font need to be measured in compara- 
ble radius, not different for each galaxy, contrary to the Einstein 
radius. In order to conclude on the favorite ongoing processes 
on ETGs since z~l as a function of their age and stellar mass, 
a joint analysis of the three lens galaxy samples would certainly 
be a first step towards a better understanding. 



5.4. Some words on the lens models 

Using lens modeling we measure the Einstein radii and related 
total masses for the lens galaxies. Doing so, we were able to 
measure their DM fractions within their Einstein radius. In ad- 
dition, the lens models were used to test the calculation of the 
external shear due to the groups around the lens galaxies. We ob- 
serve that, in every lens model, the difference between the best 
fit external shear and the shear due to the groups points towards 
the closest galaxy to the lens. We calculate that, those secondary 
galaxies need to have realistic velocity dispersions for the lens 
models to provide good fits. 

Moreover, we have measured that a catalog of high mass 
groups and clusters modeled by TIS is giving very similar total 
shear strengths and orientations than a combination of high mass 
and low mass group catalogs. Therefore we conclude that a mea- 
surement of the global external shear affecting the lens potential 
at the lens galaxy location could be made if one could combine 
: 1) a measurement of the redshift and velocity dispersion of the 
lens galaxy closest neighbor(s) plus 2) the locations, masses and 
radii of the most massive groups and clusters (such as the one 
provided by the X-ray observations in the COSMOS field) in a 
>5' radius around the lens galaxy. By fixing the external shear 
one would break an important source of degeneracy in the lens 
models. 



5.5. Conclusions 

On one hand, we have measured that the environment of lens 
galaxies is similar to that of non ETGs over a wide redshift 
range: between 0.068 and 1 . 1 3 if we put together results from the 
SLACS sample (T09, Auger 2008) and from the COSMOS sam- 
ple (this paper). On the other hand, we have built up an ensem- 
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ble of clues suggesting that the mass properties of lens galaxies 
evolve with redshift. 

Indeed, we notice that at high redshift, lens galaxies have a 
large stellar mass and a total dark matter fraction < f DM >~0.7 
(for z>0.8). On the contrary, at low redshifts, lens galaxies have 
lower stellar mass and their DM fraction are < /dm > >0.8. 

This results advocates in favor of high stellar density of high 
redshift ETGs in comparison to low redshift ETGs as suggested 
by the results of Bezanson et al. (2009). Or it could be that the 
difference between low and high redshift lens galaxy population 
is a consequence of the stellar population aging and different 
IMFs at different ages (see Napolitano et al. 2010). It could also 
be that the effects measured here are related to the lensing effi- 
ciency, which in this case would be a complex combination of 
(1) the lens population number density, (2) the source popula- 
tion distribution in space, redshift and luminosity, (3) the survey 
properties (sensitivity, band, size, angular resolution), (4) biases 
in the lens sample selection. 

To disentangle between an evolutionary or a pure lensing 
origin of the effects discovered in this study, the evolution of 
lens galaxy properties with redshift needs to be studied in ded- 
icated numerical simulations (e.g. van de Ven et al. (2009), 
Mandelbaum et al. (2009)). Whether the effects are intrinsic to 
the massive early type galaxy population or to the lensing ef- 
ficiency, they must be fully understood if one wants to study 
properly the galaxy properties gathered from lens galaxy popu- 
lations. 
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C. Faure et al.: Environment and mass properties of COSMOS lens galaxies. 

COSMOS 5921+0638 

Galaxy 2 




Fig. 8. Suite: Mass models on top of ACS images for (top to bottom) COSMOS 5921+0638, COSMOS 0038+4133 and 
COSMOS 0124+5121. The color code for the labels is given in Fig.0 North is to the top and East to the left. 
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C. Faure et al.: Environment and mass properties of COSMOS lens galaxies. 




Fig. 9. Suite: Mass models on top of ACS images for (top to bottom) COSMOS 0047+5023, J100140.12+023427.7 and 
COSMOS 5941+3628. The color code for the labels is given in Fig.H North is to the top and East to the left. 
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C. Faure et al.: Environment and mass properties of COSMOS lens galaxies. 




Fig. 10. Suite: Mass models on top of ACS images for (top to bottom) COSMOS 0050+490, COSMOS 0018+3845 and 
COSMOS 5914+1219. The color code for the labels is given in Fig.0 North is to the top and East to the left. 
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C. Faure et al.: Environment and mass properties of COSMOS lens galaxies. 

J095930.93+023427.7 




Fig. 11. Suite: Mass models for the J095930. 93+023427. 7. Top left panel: the navy blue circle located North of the image surrounds 
the group central galaxy. Top right panel: zoom on the lens galaxy and on the images. Bottom panel: zoom on the caustic (yellow 
curves) and on the favorite source position (green ellipses). The color code for the labels is given in Fig. [71 North is to the top and 
East to the left. 
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